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Abstract 

New results from the 2010 RHIC low energy program show a substantial suppression of neutral pions in central 
Au+Au collisions at both ysNN =39 and 62.4 GeV c.m.s. energies. At high pj the 62.4 GeV and 200 GeV data 
follow the same suppression pattern. On the other hand, otherwise successful pQCD predictions do not describe 
the 39 GeV data. These observations indicate that initial state effects may play a dominant role at smaller c.m.s. 
energies and at lower pj. The azimuthal dependence of the nuclear modification factor R\\ is strongly correlated 
with the (approximately elliptical) geometry of the overlap region. The dependence of R AA on the reaction plane, 
determined up to pj=20 GeV/c from 2007 high-luminosity yl'NN =200 GeV Au+Au data provides great selectivity 
among theories, and favours the ASW scenario with AdS/CFT correspondence over the pQCD-based models. 
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1. Introduction 

Neutral pion production in y/s NN =130 and 200 GeV central Au+Au collisions at the Relativistic Heavy Ion 
Collider has been found to be strongly suppressed in comparison to the expectations from properly scaled p+p col- 
lisions G]|2), while data from d+Au collisions showed no suppression or enhancement [3|, indicating that hadron 
suppression is a final state effect. This observation was one of the first convincing signatures of a strongly interacting 
partonic medium created in high energy heavy ion collisions. The nuclear modification factor R A a quantifies the 
departure from the above scaling. It is generally defined as a function of transverse momentum and pseudorapidity, 

^ T ,,)= (rAB)xd2</d/3Td?/ (1) 

where o"p° is the production cross section of n° in p+p collisions, {Tab) = N^w/cr™ 1 is the nuclear overlap function 
averaged over the relevant range of impact parameters, and N co \\ is the number of binary collisions computed with c^f'. 
The Raa, however, is also subject to initial effects such as multiple soft scattering (Cronin-effect). Measurements of 
pion production in Cu+Cu collisions at yfs^m = 22 .4, 62.4 and 200 GeV center-of-mass (c.m.s.) energies have shown 
that such effects indeed play an important role: while -^nn =200 and 62.4 GeV central data in this smaller system still 
show a suppression, pion production in 22.4 GeV Cu+Cu collisions is moderately enhanced at p T > 2 Ge V/c, virtually 
independently of centrality (4). Production of 7r°s has been measured in y.?NN =62.4 GeV and 39 GeV Au+Au 
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collisions in 2010, as part of the RHIC low energy program [5|. These data, which reach up to pj = 10 GeV/c, allow 
for the study of the onset and evolution of the suppression, and restricts theoretical models. For a better understanding 
on the interplay of the different effects of nuclear modification, other observables, such as the fractional momentum 
shift S i oss , are also measured. Azimuthal asymmetries with respect to the reaction plain carry additional information 
on the collision geometry. Recent PHENIX measurements at a/^nn = 200 GeV of n° azimuth-dependent nuclear 
modification factor Raa(Pt, A0) up to pi — 20 GeV/c [6 1 give additional constraints on energy loss mechanisms. 

2. Experiment and Analysis 

A detailed description of the PHENIX detector is given in [7 1. In the analysis of RHIC year 2007 run ys NN = 
200 GeV Au+Au data, 3.8 x 10 9 minimum bias events were used. Year 2010 run y\?NN = 39 GeV and 62.4 GeV 
analyses used 3.5 x 10 8 and 7 x 10 8 events, respectively. The Beam-Beam counters were used to determine centrality. 
The number of participants (N pm -t) and the number of binary collisions (N co n) were calculated using the Glauber 
model. In both analyses neutral pions are observed through the n° — > yy decay channel. Photons are identified with 
the Electromagnetic Calorimeter (EMCal, I77I < 0.35) 0. Only the Lead Scintillator part (PbSc, 3/47T in azimuth) 
was used in year 2010 analysis, while year 2007 high-statistics analysis used the Lead Glass (PbGl, 1 1 An) in addition. 
Extraction of the 7r° yield is detailed in (5]|6). The spectra are well described by a power-law function f(pr) = Ap T " 
at pr > 4 GeV/c, with the fitted slope parameters being n 39 = 13.4 + 0.08, n 62 = 10.4 + 0.03 and n 2 oo = 8.06 + 0.01 
for y^NN = 39, 62.4 and 200 GeV minimum bias Au+Au collisions respectively. These values are comparable to the 
corresponding values from p+p, = 13.59 ± 0.21, n'£ = 9.82 ± 0.18 and «f£ Q = 8.22 + 0.09. 

Year 2007 y^NN = 200 GeV Raa(Pt) was calculated point-by-point using p+p data from year 2005 run. At 
y^NN = 62.4 GeV/c, tt° cross section in p+p was measured in year 2006 run up to pj < 7 GeV/c. The reference was 
extrapolated into the range 7 < pj < 10 GeV/c, based on a power-law fit on the 4.5 < pj < 7 GeV/c p+p data. For 
y^NN =39 GeV, p+p data from Fermilab E706 were used, with correction for the different rapidity ranges using 
PYTHIA. Both methods import a 20% uncertainty at high p T . Decay photons from high p T pions have a small opening 
angle and, given the granularity of the EMCal, cannot be reconstructed separately. This "cluster merging" effect is 
negligible below pj < 12 GeV, but ~50% of n s are lost at pj =16 GeV, yielding a 28% error. Other main sources of 
systematic uncertainties on the pion production are PID efficiency (3-5%), yield extraction (3-5%), conversion (4%), 
acceptance (1-2%) and off-vertex n°s (1.5%). The reaction plane was determined by a combination of Muon Piston 
Calorimeter (MPC, 3.1 < \rj\ < 3.9) (TO) and the inner ring of the RxNP plastic scintillator (1.5 < \t]\ < 2.8) [11 J. The 
7r°s were classified into six bins depending on the emission angle with respect to the event plane, in order to measure 
Raa(Pt, A0). An unfolding correction was applied due to finite event plane resolution. 

3. Evolution of /? A a and Si oss with collision energy and centrality 

The 7i° nuclear modification factors /?aa(PtX measured by PHENIX in y% N — 39, 62.4 and 200 GeV Au+Au 
reactions (4][5j are shown in Fig. [T] Data from the 10% most central collisions (left panel) exhibit a strong sup- 
pression at all three energies. The mid-peripheral data (40%-60% centrality, right panel) at 62.4 and 200 GeV is 
suppressed in the whole pj range, while at 39 GeV suppression is present at low pi values, and rises slightly, though 
not significantly, above unity from pj > 3 GeV. In Fig[T] pQCD-based calculations of [12] are shown as continuous 
lines. Although the model successfully describes the 200 GeV Au+Au and Cu+Cu data, it fails at 39 and 62.4 GeV. 
A calculation with reduced Cronin-effect and with the energy loss varied by +10%, represented by shaded bands, is 
still unable to reproduce the trend in pj observed in 39 and 62.4 GeV data. The jt° Raa as a function of A^ part is shown 
in the ranges of 4 < pj < 6 GeV/c and pj > 6 GeV/c in Fig. [2] Suppression in the higher pj range is the same 
for = 200 GeV and 62.4 GeV collisions, while the 39 GeV data departs from the trend. Note that there is a 

remarkable N plLrt dependence at all energies including 39 GeV in both pj ranges (also see Fig.[TJi, and that Raa$ for all 
three energies converge at high p T . These observations, along with the radical changes of the slopes of the p T -spectra 
with c.m.s. energies and the scaling violation of the jcx-scaling effective exponent [5|, support the conclusion that 
jet quenching will be masked up to higher px, and hard scattering as a source of particles at a given pi becomes 
completely dominant only at higher transverse momentum. 

The p+p and Au+Au spectra have similar power law tails that can be fitted simultaneously with the function 
form /(/?t) = Mpi + SpiY", where Spj is the horizontal shift between the two spectra. Fractional momentum 
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Figure 1: Saa of n° in most central (0-10%, left) and mid-peripheral (40-60%, right) Au+Au collisions. Error bars are the quadratic sum of 
statistical and px-correlated systematic uncertainties, including systematics from the p+p reference. Boxes around 1 are the quadratic sum of 
overall normalization and N co n uncertainties. The continuous lines and the bands are pQCD calculations as described in the text. 
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Figure 2: Nuclear modification factor averaged for 4 < px < 6 GeV (upper 
panel) and pj > 6 GeV (lower panel). Error bars are statistical, boxes are the 
sum of px-uncorrelated and N co [\ errors. Boxes around 1 are of px-correlated 
and overall normalization errors and errors from the p+p reference. 



Figure 3 : Fractional momentum shift S i oss between Au+Au and 
scaled p+p collisions for central collisions at c.m.s energies of 
y*NN =39, 62.4 and 200 GeV as a function of the Au+Au px. 



shift Si oss = 6pj/pj is assumed to correspond to the energy loss of partons within the medium. 5i oss is shown for 
central events at three different c.m.s. energies in Fig. [3] Although the naive picture that a large suppression in R^a 
corresponds to a large S\ oss globally holds, it is to be noted that, contrary to the Raa, Si oss in 39 and 62.4 GeV is the 
same within errors below pi < 6 GeV, while in 200 GeV there is a higher momentum shift in the whole range. 

4. Collision geometry dependence of 7T° production 

There are several theoretical models for the energy loss of the high-/? T particles within the medium, such as 
the Arnold- Yaffe formalism (AMY) [14|, the Higher Twist approach (HT) lfT31 . the pQCD-based Armesto-Salgado- 
Wiedemann approach (ASW) [ 16 1, or the ASW using the AdS/CFT correspondance [ 17 1. Since the azimuthal asym- 
metry of the overlap region develops into an anisotropy of the created medium, the average pathlength a high-/?x 
hadron traverses in the medium will be correlated to the azimuthal angle A<f> between the path of the hadron and the 
reaction plane. The 2007 ysNN =200 GeV Au+Au data on azimuthal anisotropy of n production extends the pre- 
vious analysis of data from 2004 run ITPJI up to p T =20 GeV/c, and provides a stronger discrimination power on the 
models. On Fig.|4]the nuclear modification factor for the in-plane (0° < A(f> < 15°) and out-of-plane (75° < A(f> < 90°) 
7T°s are plotted and compared to the ASW scenario based on pQCD as well as the AdS/CFT correspondance. While 
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Figure 4: The data points are Raa(A0) in 20-30% centrality as a function of p~[ for in-plane (dots, upper points) and out-of-plane n° (squares, lower 
points), compared to pQCD-based ASW 1 16 (a) and ASW using AdS/CFT correspondence 1171 (b). The dotted and solid lines are the in-plane 
and out-of-plane predictions, respectively. The shaded bands around 1 are the systematic uncertainty of the normalizing ^-integrated Kaa- Boxes 
on the right side of each panel around 1 show global uncertainties. 



the out-of-plane Raa(.Pt, A</> ~ n ) lS acceptably described by all models, measured in-plane Raa(Pt, A</> « 0) clearly 
favours the ASW scenario with AdS/CFT correspondence over the pQCD-based models. The pQCD-based AMY lfT4l 
and HT lfT31l are also not supported by the data |6) . 



5. Summary 

New results from the 2010 RHIC low energy program show a substantial suppression of neutral pions in central 
Au+Au collisions at both ^Js^ =39 and 62.4 GeV c.m.s. energies. At p T > 6 GeV/c the 62.4 GeV and 200 GeV 
data follow the same suppression pattern. On the other hand, otherwise succesful pQCD predictions do not describe 
the 39 GeV data. The fractional momentum shifts 5i oss in 62.4 or 39 GeV collisions are similar below pj < 6 GeV/c, 
while S i oss is substantially higher for 200 GeV than 62.4 or 39 GeV collisions, implying a higher fractional energy 
loss. These observations indicate that initial state effects may play a dominant role at smaller c.m.s. energies and at 
lower p T , that are able to mask the effect of suppression, while the hard processes become the dominant sources of 
hadrons above pj « 6 GeV/c. 

The azimuthal dependence of the nuclear modification factor Raa is strongly correlated with the (approximately 
elliptical) geometry of the overlap region. The dependence of Raa on the reaction plane, determined up to pi-2Q 
GeV/c from 2007 high-luminosity V^nn =200 GeV Au+Au data provides great selectivity among theories, and 
favours the ASW scenario with AdS/CFT correspondence over the pQCD-based models. 
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